Magnetotransport in 2D electron systems with a Rashba spin-orbit interaction 
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The beating pattern of Shubnikov-de Haas oscillations in 2D electron system in the presence of a 
Rashba zero-field spin splitting is reproduced. It is shown, taking into account the Zeeman splitting, 
that the explicit formulae for the node position well describes the experimental data. The spin-orbit 
interaction strength obtained is found to be magnetic field independent in an agreement with the 
basic assumptions of the Rashba model. 

PACS numbers: 73.20.At,71.43.Qt,72.25.Dc, 73.61.-r 



I 

O 

o 



> 
O 



O 

in 
o 



I 

O 

o 



There has been growinff interest in the zero-magnetic- 
field spin splittingP, IJ 01 2D electron gas 
(2DEG), associated with the spin-orbit interaction 
(SOI) caused by the structural inversion asymmetry in 
heterostructures0. Apphcation of a gate voltage Hi 
is known to be the most effective method to control 
the SOI strength. These 2D systems have been sug- 
gested for application in future spintronics devices, such 
as spin-based field-effect transistors jij, spin-interference 
devices and nonmagnetic spin filters based on a 
resonant tunneling structure j^. Usually, the beating- 
pattern analysis of Shubnikov-de Haas oscillations ( 
SdHO ) jj] and the weak antilocalization method [l^ are 
used to determine the SOI strength in 2D systems. How- 
ever, the former approach is known to lead to a certain 
controversy in determining the zero-field spin splitting A. 
Namely, the spin splitting deduced from the SdHO beat- 
ing node position at finite fields 0] is different from that 
A expected for B = 0. In the present paper, this dis- 
crepancy is attributed to the contribution of the nonzero 
Zeeman spin splitting at finite fields. We support our 
idea by a rigorous analysis of the SdHO beating pattern 
caused by SOI spin splitting. The beating node posi- 
tions reported in |3| agree well with those predicted by 
the theory. Then, we demonstrate that the SOI strength 
is independent of the magnetic field. 

Let us consider a 2DEG in the x-y plane, subjected to 
a magnetic field. In the Landau gauge, the one-electron 
Hamiltonian including the Rashba spin-orbit term is 
given by 



H = 



(P 



[o-(p + eA)]n - 



(aB) (1) 



2m h'' ' ' "" ' 2 

where p is the 2D momentum; m, the effective mass; 
g, the Zeeman factor; fiB, the Bohr magneton; and, n, 
the unit vector in the z-direction. Then, a is the Pauli 
spin matrix; B, the total magnetic field; and, a, the SOI 
strength. 

It has been shown that the solution to Eq. has an 
explicit form in the case of a perpendicular magnetic field 
'B = Bz = B. The spectrum for dimensionless energy 
e — E/fi (/i is the Fermi energy) is given by 



where 77 = hcoc/fj- is the dimensionless magnetic field; 

— sR^ the cyclotron frequency; /? = ^(1 — x)j the 
term containing the Zeeman spin splitting; x = 2^ > 
spin susceptibility; and n, an integer similar to that in the 
conventional description of the Landau levels(LL). Then, 

according to Ref.0| 7 = = where S is the 

dimensionless SOI strength parameter; A = 2akp, the 
zero-field spin-orbit splitting at the Fermi energy; and 
hkp, the Fermi momentum. Usually, the typical Fermi 
energy fi ^ 80meV exceeds the SOI- induced splitting 
A ~ ImeV ( see Q ), and, therefore (5^1. It is notewor- 
thy that the conventional spin-up(down) energy states 
associated with n-th LL number correspond to and 
states respectively. In the absence of SOI, Eq.Q 
reproduces well-known LL energy spectrum. 

In contrast to the conventional formalism extensively 
used to find the low-B rnagnetoresistivity, we use the al- 
ternative approach 0, ITsl which allows to resolve 
magnetotransport problem in both the SdHO and Inte- 
ger Quantum Hall Effect (IQHE) modes. Moreover, this 
method was successfully used in a recent paper [13 to 
reproduce the SdHO beating structure in the presence of 
the zero-field valley splitting (Si-MOSFET 2D system), 
and in both the crossed- and tilted- field configurations. 
Following the argumentation put forward in Ref . [l^ , well 
above the classically strong magnetic field range WcT ^ 1, 
where r is the momentum relaxation time, 2DEG can be 
assumed dissipationless in strong quantum limit when the 
cyclotron energy fiiOc exceeds both the thermal energy 
kT and the energy related to LL- width h/Tq. Here, Tq is 
the quantum relaxation time. Under the above assump- 
tions a XX 1 Pxjr. ^ 0. Nevertheless, routine dc measure- 
ments yield [ig the finite magnetoresistivity associated 
with a combination of the Peltier and Seebeck thermo- 
electric effects. Within the scenario suggested 0, we 
obtain the above magnetoresistivity in the form 



p = Pyx- 



'■2D 



(3) 



(2) 



where a2D is the 2DEG thermoelectric power; — 
Nec/B, the Hall resistivity; TV = -(|^)^, the 2D den- 
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sity, n = -fcTrX;in(l + exp(ii^)), the thermody- 



namic potential; F 



= the zero-width LL density of 
the Lorentz number; ks, the Boltz- 



mann constant. In fact, the 2D thermoelectric power in 
strong magnetic fields is a universal quantity [l^, pro- 
portional to the entropy per electron: a2D — where 
S = -(fp)^ is the entropy. Both S,N, and, there- 
fore, a2D,P are universal functions of the dimensionless 
temperature £, — ^ and the magnetic field rj — 2/^, 
where v = Nq/T is the conventional filling factor, and 
-^0 = ^Fp'M is the zero- field density of the strongly degen- 
erate 2DEG in the absence of a SOI-induced splitting. 

Using the Lifshitz-Kosevich formalism and, then, ne- 
glecting finite LL-width( h/xg — > ), we derive in 
Appendix asymptotic formulae for Q, and, hence, for 
N, S, Pyx , p, which are valid at low temperatures and 
weak magnetic fields ^,r] <^ 1: 

N = m^Fo{m) + 2n^No±'-^^^Riv), (4) 

OO 

S^So- in^^ksNo cos(27rfc/77)i?(r,), 

fc=i 

where So = kBNo{2^Fi{l/() ~ Fo{l/()) is the entropy at 
5 = 0; F„(z), the Fermi integral; and $(z) = ^~ sinhl'i)'^ ■ 
At -B = both the thermopower and 2D density are 
constants, i.e. a2D = a""?"' ^ — ^o, hence the magne- 
toresistivity is given by zero-field asymptote p — ^ g . 
According to Eq.Q, for actual first-harmonic case( /c = 1 
) the magnetoresistivity can be viewed as the zero-field 
background, on which the rapid SdHO modulated by 
long-period beatings( see Fig[51) are superimposed. It's 
worthwhile to mention that at the beat nodes( i.e. when 
the form-factor at fc = 1 vanishes ) the magnetoresistiv- 
ity is given by zero-field asymptote. This is not, however, 
the case of low temperatures and(or) high magnetic fields 
when the high-order terms(fc > 1) in Eq.l^J may de- 
termine the amplitude of magnetoresistivity at the beat 
nodes. It turns out that the data reported in point to 
the above feature. 

We now analyze in detail the form-factor R{rj){ see 
Appendix ) which determines the beating pattern of S, N 
and, hence, p. For the actual first-harmonic case (i.e., 
A; = 1), the beating nodes can be observed when R{rj) = 
or 



(5) 



where we neglect the small quadratic term S^/irj^ ^ 
d/r]"^ evaluating Eq.©. Then, j = 1,3.. is the beat- 
ing node index. We emphasize that the first node cannot 
be observed in experiments, performed, for example, in 
Ref. Indeed, for real 2D Ina;Gai_^As/Ino.52Alo.48As 
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FIG. 1: Zero-field SOI splitting at the Fermi energy vs node 
index j, deduced from the experimental data Q , with the help 
of the node condition specified by Eq.llKJ. Dotted lines( from 
top to bottom ) represent the mean values Aq for samples 
B,A,C respectively. 

TABLE I: Transport data( at 4.2K ) and the Zero-field spin- 
orbit splitting at Fermi energy for Ina,Gai_a: As/Ino.52Alo.48As 
2D system reported in Ref.Q 



Sample(x) 


/io X 10*,cmVVs 


n X lO^^cm"^ 


pi,meV Ao,meV 


A(0.65) 


13.4 


1.75 


78 2.34 


B(0.60) 


9.5 


1.65 


74 2.57 


C(0.53) 


6.8 


1.46 


65 1.63 



system ( m — 0.049mo, g ~ 4 ) we find (3 — 0.45, and, 
therefore Eq.© cannot be satisfied for j = 1. With the 
help of Eq.©, we analyze the nodes, reported in for 
three different samples, and then plot the dependence of 
the zero-field SOI splitting at the Fermi energy A against 
the node index( see Fig^), starting from j = 3. For these 
samples A is nearly constant within the actual range of 
the magnetic fields, therefore we obtain the respective 
mean values Aq denoted in Tabled Note that the minor 
deviation of A with respect to its mean value in high-field 
limit ( low- index nodes ) can be associated with possible 
magnetic field dependence of the g-factor. In contrast, 
the non-parabolicity effects [l^ seem to be irrelevant 0| 
for the actual low-field case B < IT. 

We emphasize that the node condition similar to Eq. Q 
was previously discussed in literature. Following the 
analysis done in Ref.0, the nodes occur when the spin- 
orbit-split subbands are shifted one with respect an- 
other by half a period at the Fermi energy. Namely, 
1 =^ ei 



i+s+i)/^! where s = 0,1,2... cor- 
responds to the node index as j = 1 -I- 2s. For actual 
high LL-number case n ^ 1 this condition reproduces 
Eq.©. 

Let us discuss the conventional method often used 
to extract the zero-field SOI splitting at the Fermi en- 
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FIG. 2: SdHO beating pattern calculated with the help of 
Eq.lEll at fc = 1 for sample A(a; = 0.65) ;3j: = 1.75 * 10^^ 
cm"^ m = 0.049mo, p = 4, Aq = 2.34meV, T = 1.6K. Ar- 
rows show the beating nodes at j = 3,5... Zero-field asymp- 
tote is represented by dotted line 

ergy. According to phenomenological arguments put for- 
ward by Das et al [H, the nodes may occur when 

cos (^T^^;^^ = or Afot ~ i^fiuJc, where the to- 
tal spin sphtting at the Fermi energy between spin- 
down e^^i and spin- up e+ states yields Atot = fn^c — 

{2(3hwcY + A^. As expected, the total spin split- 
ting Atot coincides with the zero- field —A and the Zee- 
man spin splitting in low( high ) magnetic field 
limit respectively. With the help of the dimensionless 
units the node condition suggested by Das et al reads 

P"^ + ^ = ^^2'^ 7 hence, reproduces our result if one 
selects set at j > 1. We argue that straightforward 
procedure ( see Fig^ ) used to extract Aq is, however, 
preferable compare to zero-field extrapolation method 
suggested in Ref.|3,0|. Indeed, for low-density samples 
and(or) under the temperature enhanced conditions the 
SdHO amplitude is suppressed, hence, the low-field nodes 
become hidden. In this case the zero-field extrapolation 
method Q may lead to a subsequent errors. 

Let us now reproduce( see Fig|21 ) the SdHO 
beating pattern with the nodes occurred in a typi- 
cal sample( sample A(a; — 0.65) in 3]) at i? — 
0.873; 0.46; 0.291; 0.227; 0.183; 0.153T using Eq.^, and 
previously extracted value of zero-field SOI splitting 
Ao = 2.34meV. It's worthwhile to mention that our 
results differ with respect to those, which can be ob- 
tained within the conventional formalism in the following: 
(i) the low-field quantum interference, classical magne- 
toresistivity and 3D substrate parallel resistivity's] back- 
ground are excluded within our approach; (ii) in contrast 
to conventional SdHO analysis, our method determines 
the absolute value of magnetoresistivity, and, moreover, 
can lead to a gradual transition [Tell from the SdHO to 



the IQHE mode. 

We argue that the noticeable increase in SdHO am- 
plitude was observed at i? ~ 0.3T. This value satis- 
fies the criterion of the classically strong magnetic field 
since lOcT = 4 while the corresponding cyclotron energy 
hujc = 8.2K correlates with that ~ 9.8K expected from T- 
dependent SdHO-damping factor, i.e. when 27r^^/77 ^ 1. 
We conclude that the energy associated with LL width 
~ h/Tq is less or at least equal to the thermal energy. 
The above estimates point to validity of zero-width LL 
model in this particular case. Nevertheless, since both 
the temperature and finite LL width known to suppress 
the SdHO amplitude in a rather similar manner, we es- 
teem reasonable to reproduce in Fig|21the SdHO beating 
pattern using somewhat higher temperature T = 1.6K 
than that T = 0.5K reported in 0. 

Note that our approach provides a correct number of 
oscillations between the adjacent nodes. For example, 
the number of oscillations confined between j = 3,5 
nodes (37) correlates with that (35) observed in 0. A 
minor point is that our approach predicts a somewhat 
lower amplitude of SdHO, compared with that in the ex- 
periment 1^. For example, for j = 3 node( B — 0.873T in 
Ref.0 ) we obtain p — 0.0035Ohm. Actually, one would 
expect the same order of magnitude for SdHO amplitude 
between the proximate nodes ( see j = 3, 5 in FigI21). 
Our estimation is, however, less than both the absolute 
magnetoresistivity 40Ohm at B = 0.873T and SdHO am- 
plitude ~ 50hm reported in Ref.P]. 

In conclusion, we demonstrated the relevance of the 
approachjlij regarding the beating pattern of SdHO 
caused by Rashba spin-orbit interactions. Taking into 
account the Zeeman splitting, the rigorous analysis of 
experimental data 's'l suggests a B-independent strength 
of the Rashba SOI. The above findingis consistent with 
the general theoretical assumptions j^. Our approach 
can be helpful for estimation of the SOI strength. 
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ported by the Russian Foundation for Basic Research 
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APPENDIX 

Using the conventional Poisson formulae 

00 °9 00 °9 

J2'p{n)= / (p{n)dn + 2ReY^ / (p{n)e^"''"dn, (6) 

where tuq — 1 < a < mg , mg the initial value of the sum- 
mation, the thermodynamic potential can be represented 
as the sum 51 = Oq + of the zero-field and oscillating 
parts as follows 

no = -No^^eFlil/0. (7) 
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In 1 + e" 



where Fn{z) is the Fermi integral. For simphcity, we omit 
the SOI-induced sphtting in zero-field term flo because 
5^1. The special interest of the present paper is in the 
oscillating term fl^ of thermodynamic potential, which 
can be strongly affected by spin-orbit-split subbands(±). 
After a simple integration by parts, the oscillating term 
yields 



27rfc 



J2-Kikn 



1 + e « 



-de 



(8) 



Using Eq.lO, for a certain energy we calculate the ac- 
tual high-order LL-likc numbers, associated with both 
the spin-orbit-split subbands as 



7^2 



V 



(9) 



It should be noted that the integrand equation in Eq.© 
is a rapidly oscillating function, which is, in addition, 
strongly damped when e > 1. The major part of the 
magnitude of the integral results from the energy range 
close to the Fermi energy, when e 1. Therefore, n^(e) 
can be regarded as smooth functions of energy, and, 
hence, can be re-written as n 



± _ 
± 



where we use the designation Ui 



"(1) . Under the 
above assumption, we can change the lower limit of in- 
tegration to — cxD and then use the textbook expression 

/ f-f^dy = sinh(lk) the integral of the above type. 

— oo 

Finally, the thermodynamic potential yields 



k=l 



tdn^ 



cos(7rfc(n^ -I- ))R{'n) 
Tk sinh Tk 



(10) 



, ~ i is valid for the ac- 
1 »; 

=*= ^ 1, and 



where we assume that 

tual case of high- order Landau levels ^ 
Tk = 27r^^A:/7/ is a dimensionless parameter related to T- 
damping of SdH amplitude. Then, i?(r/) = cos(7rfc(nf — 
is the form-factor. The oscillatory part of the 
thermodynamic potential consists of rapid oscillations 
cos(7rfc(n^ -|- cos(27rfc/77), on which long-period 



beatings governed by the form-factor are superimposed. 
As expected, the form-factor is reduced in absence of SOI 
to a field-independent constant R{ri) = cos(27rfc/3), and, 
therefore, the beating structure is absent. Using the con- 
ventional thermodynamic definition, we can easily obtain 
both the entropy and the density of 2D electrons, speci- 
fied by Eq.Q). 
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